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I. INTRODUCTION 

An ever increasing amount of evidence indicates the existence of large amounts of dark 
matter in the Universe Jl|. Despite this overwhelming evidence, the exact nature of the 
dark matter remains a mystery. Numerous candidates have been proposed, including both 
baryonic and nonbaryonic matter || . Observations reveal that some of the dark matter in 
the Galactic halo is baryonic, consisting of MACHO's ||; however, present data indicates 
that MACHO's cannot account for all of the dark matter implied by the Galactic rotation 
curve 0. Furthermore, a number of arguments based upon large scale motions in the 
Universe and large scale structure formation indicate that f]«l, which is far in excess of 
the bounds on ^naryon^ 2 < 0.026 arising from cosmic nucleosynthesis ||. All considerations 
point toward non-baryonic matter comprising a sizable fraction of the Universal density. If 
this is true, what is the dark matter? 

Among the best motivated, and hence highly favored, of the nonbaryonic dark matter 
candidates is the lightest supersymmetric particle (LSP). Experimental and theoretical con- 
siderations indicate that the LSP is a neutralino, x, consisting of a linear combination of 
the supersymmetric partners of the photon (7), the Z (Z), and 2 Higgs bosons (Hi and H2). 
Note that the 7 and Z are themselves linear combinations of the supersymmetric partners of 
the neutral W (IU3) and B (B) bosons, hence the neutralino composition is typically written 
as 

X = Z t B + Z 2 W 3 + Z Z E X + Z A H 2 . (1) 

The motivation for supersymmetry (SUSY) arises naturally in modern theories of particle 
physics |l|||, although the x's potential as a dark matter candidate was not realized until 
later [7|]. For a very large region of SUSY parameter space, neutralinos provide densities 
that are in accord with the measured value of Q, thus explaining the dark matter. The x is 
also detectable in at least two ways: indirectly, through the products of XX annihilation in 
the Sun, Earth, and Galactic halo [|l]f|, or directly, via elastic (and inelastic |§) neutralino- 
nucleus (xN) scattering in a detector p|,p!0|. In either case, the elastic xN scattering cross 
section is an essential ingredient. In this paper we discuss nuclear structure calculations 
relevant to xiV scattering for several nuclei which are primary constituents of many current 
and planned direct detection experiments [p|,pT|-[T5|] . 

Physics at three distinct energy scales governs xN scattering. The composition and mass 
of the x, and hence its interaction with quarks are fixed near the electroweak scale. The 
interaction of neutralinos with protons and neutrons is determined by the quark distribution 
(both spin and density) within the nucleon, which is determined at the QCD scale. At 
the modest momentum transfers available to dark matter neutralinos the x interacts with 
the entire nucleus, not individual nucleons within it. Thus, nuclear structure plays an 
important role in determining the xiV cross section. The uncertainties in the electroweak 
scale physics (the SUSY part of the problem) are typically handled by considering large 
sweeps through SUSY parameter space 0,0- The QCD scale physics is currently the focus 
of much study and the relevant nucleon matrix elements continue to be measured with 
high precision [|TJ]]. The necessary nuclear physics is not measurable for most nuclei but is 
amenable to calculation through a variety of methods. Here we apply the nuclear shell model 
to the nuclei I, 129 Xe, 131 Xe, 125 Te, and 23 Na (this last nucleus we discuss in an appendix) 



in order to provide a consistent and correct set of nuclear input physics for determining the 
xN cross section. 

The \N scattering cross section has two distinct terms: a spin- independent, or scalar, 
term, and a spin-dependent, or axial term. It has been well established that for nuclei with 
A > 30-50 (A = number of nucleons), the scalar piece of the interaction tends to dominate 
the xiV scattering rate; however, there are significant regions of parameter space where this 
is not so and the axial rate dominates |j]|R| . The importance of understanding the axial xiV 



interaction is amply demonstrated by a recent SUSY interpretation of a Fermilab scattering 
event |lij|| . These papers claim that the \ might be an almost pure higgsino whose couplings 
to ordinary matter are completely dominated by the axial part. In this paper, we shall deal 
with the axial \N interaction. The relevant nuclear physics for the scalar interaction is 



well approximated by a fairly simple form factor, suitable for all nuclei p9| , p0| . The axial 
response is far more complicated and requires detailed nuclear models. 

II. NUCLEAR PHYSICS INPUT 

A variety of nuclear models have been used to calculate the axial response of nuclei 
used as targets in dark matter detectors. The conventional nuclear shell model J2T| has 
proven highly successful at accurately representing this response when a reasonable nuclear 
Hamiltonian is used in a sufficiently large model space |]22|-p5|]. Until recently, both of 
these ingredients have been absent for nuclei in the 3s2<ilg7/2l/iii/2 shell, including most 
of those included in this study. With recent advances in computer power and storage, we 
can now construct model spaces that contain most of the nuclear configurations that are 
likely to dominate the spin response of nuclei such as 12 I. Coupled with this ability to 
perform sufficiently large calculations is the recent development of several realistic nucleon- 
nucleon (nn) potentials [p6|,2"7J1. These potentials can then be converted into suitable nuclear 



interaction Hamiltonians via the G- matrix/Folded Diagram technique [p6 ]. In this paper 



we consider two such nuclear interactions, one using the Bonn A p6| and the other the 



Nijmegen II [27] nn potential. The diagonalization of the Hamiltonian was performed using 



the shell model code ANTOINE 



A. The Hamiltonian 

The residual nuclear interaction based upon the Bonn A potential has been described 
in Ref. p6| . This Hamiltonian has been derived for the model space consisting of the 
1^7/2, 2g? 5 /2, 3si/2, 2<i 3 /2, and l/in/2 orbitals (which we use in this study). It was originally 
derived to describe light Sn isotopes (A « 102-110) which have no protons in the space. In 
order to find good agreement with observables for nuclei with A « 130, the single particle 
energies (SPEs) were adjusted. We made an initial guess at the SPEs based upon the excited 
state energy spectra of nuclei with either a single neutron hole in the space ( 131 Sn) or a single 
proton in the space ( 133 Sb and 125 Sb). These initial SPEs were then used in conjunction 
with the Two Body Matrix Elements (TBMEs) of the interaction to calculate observables 
for the nucleus I. We varied the SPEs until reasonable agreement between calculation and 
experiment was found for the following 127 I observables: the magnetic moment (/i), the low 



lying excited state energy spectrum, and the quadrupole moment (Q2o)- This procedure 
is similar to that performed in Ref. Iffil - The magnetic moment is extremely important, 
as it is the observable most closely related to the \N scattering matrix element and has 
traditionally been used as a benchmark of a calculation's accuracy. In Fig. [I], we show the 
final SPEs used in our calculations. In Table | we show the final calculated values of \x and 
Q20 for 127 I vs. the experimental values; agreement is excellent. Once the SPEs are specified, 
we have a reasonable Hamiltonian to use for the nuclei we are studying. 

In order to test the sensitivity of our results to the nuclear Hamiltonian, we have also 
examined another one, derived from the Nijmegen II nn potential [27| . We have used 
the codes and methods described in []29| to convert the potential to a usable shell model 
interaction. The procedure is similar to that used for the Bonn A force. The two sets of 
TBMEs are generally similar but significant differences do exist. We initially used the same 
set of SPEs as above but found that a significant lowering of the 1(77/2 SPE was necessary in 
order to find agreement with the observables. The SPEs and comparisons with observables 
for this force are shown in Fig. [I] and Tables | and [TJ. 

B. The Model Space 

To perform a full basis, positive parity, calculation of the 127 I ground state properties 
in the space consisting of the 1(77/2, 2<i 5 / 2 , 3si/2, 2rf 3 / 2 , and l/in/2 orbitals, we would need 
to have basis states consisting of roughly 1.3 x 10 9 Slater Determinants (SDs). Current, 
state of the art, calculations (including those presented here) can diagonalize matrices with 
basis dimensions in the range 1-2 xlO 7 SDs; clearly severe truncations of the model space 
are needed. Fortunately, given the size of the model spaces that can be treated, a truncation 
scheme that includes the majority of relevant configurations can be devised. 

Our scheme is best understood by viewing Fig. |1|. As a base configuration, we have for 
protons: ( 1*77/2 2rf 5 / 2 ) 3 (i.e. a total of 3 protons spread among the 1(77/2 and 2d 5 / 2 orbitals) 
and for neutrons: (l(77/22<i5/2) 14 + (3si/ 2 2<i3/2) 6 + (ih u / 2 ) 4 . We then assign the following 
values of the "jump" to each orbital: j ump( lg 7 / 2 , 2<i 5 / 2 ) = 0, jump(3si/ 2 , 2d 3 / 2 ) = 1, and 
jump(l/in/ 2 ) = 2. The differences in these values is the cost of moving particles between the 
different (sets of) orbitals. Hence, to move a proton from the 1(77/2 to the 2d§/ 2 costs nothing 
while moving one from the 1(77/2 to the 3si/ 2 would cost 1 unit of jump (to the lh u / 2 would 
cost 2 units). It would cost 2 units of jump to move 2 neutrons from the 2d 3 / 2 to the l/in/2, 
etc.... All that remains is to specify the total amount of jump available. In our truncation, 
we allow protons up to 3 units of jump, neutrons up to 4 units, and a total of up to 4 when 
adding the jump used by the protons plus neutrons. Thus, if the protons remain in the 
1(77/2 and 2d 5 / 2 orbitals (as they tend to do), the following neutron configurations are al- 
lowed: (lg 7/2 2d 5/2 ) 14 + (3s 1/2 2d 3/2 ) 6 + (l/i 11/2 ) 4 , (lg 7/2 2d 5/2 ) u + (3s 1/2 2d 3/2 ) 4 + (l/ni/2) 6 , 
(lg 7/2 2d 5/2 ) 13 + (3s 1/2 2d 3/2 ) 5 + (l/m/2) 6 , (l27/224/ 2 ) 14 + (3si/22d 3/2 ) 2 + (lh u/2 ) 8 , and 
{lg7/22d 5 / 2 ) 12 + (3s 1 / 2 2d 3 / 2 ) 6 + (l/in/2) 6 - If 1 or 2 protons are excited out of the 1(77/2 
and 2g?5/2 orbitals, the last 2 neutron configurations are not allowed. In this truncation, the 
m-scheme dimension of the 127 I model space is about 3 million SDs. 

Our results indicate that this space is more than adequate to describe the ground state 
properties of the nuclei considered. As mentioned above, our calculation of the observables 
agrees well with experiment. The major potential problem with this model space would be 



if it failed to allow enough neutron excitations into the l/in/2 orbital . It allows at most 8 
neutrons out of a possible 12 in that orbital. In Table [TlT| we present the occupation numbers 
for 127 I. We see that our interactions do not seem to prefer excitation of more than one extra 
neutron pair to the l/in/2- Most configurations have six neutrons in that orbital, while eight 
are allowed. Hence, our model space is more than adequate. 

For the two Xenon isotopes considered (A = 129 and 131), we have used exactly this 
truncation scheme. For 125 Te we used this scheme and also one where the total jump and 
total neutron jump was 6 (instead of 4). Very little difference was noticed for the two 
truncations. In this paper we present the results for the larger truncation since it should be 
slightly more realistic. 

III. RESULTS 

A. The Zero Momentum Transfer Limit 

Neutralinos in the halo of our Galaxy are characterized by a mean virial velocity of, 
v ~ (v) ~ 300 km/sec = 10 _3 c. The maximum characteristic momentum transfer in %iV 
scattering is q max = 2M r v where M r is the reduced mass of the \N system. If the product 
q max R is small (<^C 1), where R is the nuclear size, the matrix element for spin-dependent 
xN scattering reduces to a very simple form |],[19| 

M = C(N\a p S p + a n S n \N)-a i (2) 

where 

Si = 5^Si(k), i = p,n (3) 

k 

is the total nuclear spin operator, k is a sum over all nucleons, and a p , a n are x-nucleon 
coupling constants which depend upon the quark spin-distribution within the nucleons and 
on the composition of the \. Much of the uncertainties arising from electroweak and QCD 
scale physics are encompassed by a p and a n . The normalization C involves the coupling 
constants, masses of the exchanged bosons and various LSP mixing parameters that have 
no effect upon the nuclear matrix element. Eq. (2) has often been written as 

M = CA(N\ 3\N) ■ s^ (4) 

with 

_ (N\a p S p + a n S n \N) (N\(a p S p + a n S n ) ■ 3\N) 

(N\3\N) J(J + 1) ' U 

Examples of the full %iV cross section can be found in Refs. |I,19j22]|. 

Equations (2-5) show that the \ couples to the spin carried by the protons and the 
neutrons. The matrix element (2) is similar to the magnetic moment operator: 

fi = (N\g s n S n + g l n L n + g s p S p + g l p L p \N). (6) 



The free particle g-factors are given by: g s n = —3.826, g l n = 0, g^ = 5.586, and g l = 1 (in 
nuclear magnetons). Given the similarities of eqs. (2) and (6), it is no surprise that // is 
often used as a benchmark on the accuracy of the calculation of S p and S n in A. We follow 
that prescription as well. In the following section we will briefly outline some problems with 
this procedure. 

In Table |T| we present the values for: S p , S n , L p , L ra , and \x that we calculate for each 
Hamiltonian for the nuclei 127 I, 129 Xe, 131 Xe, and 125 Te. We also include the experimentally 
measured magnetic moment, it is apparent that agreement is quite good for all nuclei. A 
number of other calculations of these quantities appear in the literature, and we include 
summaries of these calculations in the table as well. The following abbreviations are used 
for the various nuclear models: Bonn A = our calculation using the Bonn A derived force, 
Nijmegen II = our calculation using the Nijmegen II derived force, OGM = the Odd Group 
Model |30|], IBFM = the Interacting Boson Fermion Model [31|, ISPSM = the Independent 



Single Particle Shell Model |[32| , TFFS = the Theory of Finite Fermi Systems |p3fl , and 
QTDA = the Quasi Tamm-Dancoff Approximation ||20|| . In most previous experimental 
analyses, the OGM values have been used [fP2"UT3"| . 



Examining our results for S p and S n in Table H and comparing them to results from other 
nuclear models reveal several interesting facts. In almost every instance, our results show 
that the spin |Sj| (i — p,n) carried by the unpaired nucleon is greater than that found in 
the other nuclear models (except for the ISPSM, where |S*| is maximal). Despite our larger 
values for | S$| , our calculations have significant quenching of the magnetic moment and are 
in good agreement with experiment in all cases (see the later section on quenching). The 
reason that we find larger values of |S»| for the odd group is due to the fact that we allow 
more excitation of the even group of the nuclei; allowing them to be a major contributor 
to the total nuclear spin: J = S p + S n + L p + L n . The naive expectation for L n in the 
Bonn A calculation of 127 I is zero. We find L n = 0.779, L n is responsible for over 30% of 
iodine's total angular momentum (J = |). This explains both the large quenching of fi 
(L n does not contribute to \i since g l n = 0) and the large value of S p found. We note that 
most previous experimental analyses used the OGM value for 127 I, S p = 0.07. Our results 
give a factor of ~ 20 increase in iodine's sensitivity to spin-dependent scattering over that 
previously assumed. Due to the form factor suppression (discussed below) a sodium iodide 
detector's [JTT1 , |T3"1 spin response is still dominated by 23 Na but not to the extent previously 



thought. For the remainder of the nuclei considered Table || also reveals increased scattering 
sensitivity, although the factor of increase is much more modest. 

B. Quenching and Uncertainties 

As we noted earlier, the comparison of the computed magnetic moment vs. the ex- 
perimental value has been used as the primary (and in some cases, only) indicator of a 
calculation's reliability. This seems quite reasonable in light of the similarities between the 
matrix elements in eqs. (2) and (6). This prescription is not without several potential 
problems |2^,[34|]. Not only does \x depend upon the orbital angular momentum Lj but the 



spin angular momentum Sj is subtly different. The j(N matrix element (6) results from the 
non-relativistic reduction of the axial- vector current. Because of this, it is not strongly af- 
fected by meson exchange currents (MECs). The magnetic moment's spin operators, S», are 



a result of the non-relativistic reduction of the vector current. They can be strongly affected 
by MECs [23,34]- The effects of MECs upon \x is typically lumped together with several 



other effects to give effective ^-factors |35|,|36|]. Unfortunately, there is no hard and fast rule 
as to what effective ^-factors are the best. We have chosen to remain with the free particle 
^-factors. As an example of the potential uncertainties this ambiguity leads to, we have 
also included, in Table [HJ the calculated magnetic moment for our nuclei using a reasonable 
set of effective ^-factors. The "quenched" magnetic moments are the values in parentheses 
in the table and the effective (/-factors used are: g s n = —2.87, g l n = —0.1, g s v = 4.18, and 
g l — 1.1. The table shows that these ^-factors do little, overall, to improve the concordance 
between calculation and experiment. 

A related concern involves the quenching of the (isovector) Gamow- Teller (GT) ^-factor, 
gA p2p4fl ■ The spin term of the GT operator also comes from the axial vector current 
and thus is closely related to the spin operators in eq. (6). Its is well established that 



most nuclear model calculations of GT strength require a reduction of g^ of order 20% 
Whether this quenching of g^ should also be applied to a± (the isovector x-nucleon coupling 
constant) is unknown |53 ]. Since there is no real guidance, and our magnetic moments 



agree well with experiment, we do not believe that any extra quenching of the spin matrix 
elements (or equivalently the coupling constants a and ai) is desirable for these nuclei when 
calculating \N scattering rates. Nonetheless, as pointed out in Ref. 131] it is useful to keep 



these potential uncertainties in mind when calculating scattering rates. 



C. Finite Momentum Transfer 



When the LSP was first proposed as a viable dark matter candidate, its preferred mass 
was between 5 and 10 GeV 0. With a mass of this order and a typical galactic halo velocity 
(v ~ 10~ 3 c), the neutralino's total momentum (q ~ M r v ~ lOMeV) was small compared to 
the inverse of the nuclear size (1/R ~ 1/1 fm ~ 200 MeV) and the zero momentum transfer 
limit was appropriate for studies of xiV scattering. Since then, experiments at accelerators 
have pushed the allowed x mass, m^, to larger values (there are ways around this if some of 
the theoretical assumptions are relaxed |J7|), and it has been shown that heavy x's are just 
as viable as a dark matter candidate as the lighter ones [^,^] . As m% becomes larger than 
a few 10's of GeV the product qR starts to become non-negligible and finite momentum 
transfer must be considered for heavier nuclei. 

The formalism for elastic xN scattering at all momentum transfers has been developed in 
Refs. 



[0,0. Here, we follow precisely the definitions used in |22 
from our definitions to those used in Ref. 



It is a simple matter to go 



. The formalism is a straight forward extension 
of that developed for the study of weak and electromagnetic semi-leptonic interactions in 
nuclei ||Q]. The differential xiV cross section is given by 



da 

dq 2 



8G F 



;S(q), 



(2J+l)v 2 
where S(q) is the spin structure function 

S(q)= E (l(^||^ 5 (g)||iV)| 2 + |(iV||£i(g)||iV)| 2 ). 



(7) 



(8) 



Lodd 



T el5 (q) and C 5 (q) are the transverse electric and longitudinal multipole projections of the 



axial- vector current operator | 40| . The double vertical lines imply that these are the reduced 



matrix elements of these operators. For their explicit form in the \ context, see |J19| , |20"| , |2"2" 
In the limit of zero momentum transfer S(q) reduces to 

9 7-1-1 
S(0) = ^^A 2 J(J + 1). (9) 

7T 

The reduced matrix elements of the multipoles in eq. (8) are easily evaluated in the har- 
monic oscillator basis in the nuclear shell model EO] . With the exception of the calculation of 



the 27 A1 structure function in []23| all calculations of S(q) have used bases of these harmonic 
oscillator wave functions. In this paper, we have used the more realistic Wood-Saxon wave 
functions to evaluate eq. (8). To specify the wave functions, we use the parameters recom- 



mended in [41|. We have used the codes from [42] to calculate the actual wave functions. We 



have also calculated the Bonn A structure function for 127 I using harmonic oscillator wave 
functions. The differences in the two prescriptions are significant at very large momentum 
transfers but are minor at most relevant values of of the momentum transfer (q) . 

It is useful (and traditional) to use the isospin convention instead of the proton-neutron 
formalism when discussing \N scattering at finite momentum transfer. Writing the isoscalar 
coupling constant as a = a n + a p and the corresponding isovector coupling constant as 
a\ = a p — a n we may split S(q) into a pure isoscalar term, Soo, a pure isovector term, Sn, 
and an interference term, Sqi, in the following way: 

S(q) = alSo (q) + a 2 S n (q) + a a x S 01 (q). (10) 

Using this decomposition of S(q) it is a simple matter to derive the structure function for a 
X of arbitrary composition. 

Two factors contribute to the maximum allowed momentum transfer. As m^ becomes 
much greater than the nuclear mass, m^v, the reduced mass asymptotes to M r — > m^. 
Also, the x' s have a Maxwellian velocity distribution in the halo and some will possess 
velocities significantly greater than (v) ~ 10 _3 c. A maximum velocity of v max ~ 700 km/sec 
(slightly greater than Galactic escape velocity |14]]) implies maximum momentum transfers 
of q m ax{A ~ 127) ~ 550 MeV. This value is not small compared to the inverse nuclear size. 
In a harmonic oscillator basis, the fiducial nuclear size is set by the oscillator parameter, 
b = lfmv4 1//6 = (l/197.327MeV)yl 1 / 6 . In order to maintain contact with previous literature 
2"2|,|2"5[ we retain b as the size parameter in our Wood-Saxon evaluations of S(q). We do, 



however, use a slightly better, empirical, parameterization of b jKJ: b = (41.467/Ziu;) 1 / 2 fm 
with huj = 45A" 1 / 3 - 25A~ 2 / 3 MeV. Hence, we have values near b(A = 127) = 2.282 fm 
= 1/86.47 MeV for the nuclei in this study. We parameterize all of our structure functions 
in terms of y = (qb/2) 2 . For ?/ <C 1 the effects of finite momentum transfers are small; 
for y > 1 the effects are quite noticeable. For these nuclei y max = {o max b/2) 2 ~ 10 3> 1, 
hence nuclear form factors are extremely significant. These extremely large values of y are 
only valid for extremely massive x' s moving near escape velocity. A more realistic x, with 
m^ = 100 GeV moving at (v) would have y ma x — 0.4. 

In order to cover all of the relevant \ parameter space, we have evaluated the structure 
functions all the way to y — 10 for the nuclei studied. This presents a problem, in that it 
has become standard to present structure functions as polynomials in y of order 6 or less. 

8 



(A structure function of this sort can easily be incorporated into the code Neutdriver of Ref. 
[ij.) We could find no suitable fits of this form valid out to values of y — 10. We have 
addressed this in two ways. In appendix B we present fits of the structure functions S 00 , 
Sqi, and Sn as 6th order polynomials in y. These fits are only good for values of y < 1. 
They should NOT be used beyond this value as they give meaningless results. In order 
to accurately represent S(q) at all relevant momentum transfers, we have had to resort to 
a somewhat more complicated functional form. In a harmonic oscillator basis, the matrix 
elements of the operators T el5 (q) and C 5 (q) are precisely represented as polynomials in y 
times a factor of exp(-y). (The isovector Goldberger-Trieman term in C 5 (q) complicates this 
slightly.) Using this form as a guide we have have fit the structure functions as 8th order 
polynomials in y times a factor or exp(-2y) This form has proven adequate to accurately 
describe the structure functions for 127 I and 131 Xe. A slightly more complicated form with a 
term added to mimic the effect of the Goldberger-Trieman term was required for 129 Xe and 
125 Te. As an example, we present the fit for the term Soo(g) for the Bonn A calculation of 

127j. 

S 00 (y) = e" 2?/ (0.0983393 - 0.489096?/ + l.U02y 2 - 1.47168?/ 3 + 1.1717?/ 4 

-0.564574?/ 5 + 0.158287?/ 6 - 0.0238874?/ 7 + 0.00154252?/ 8 ) (11) 

We relegate the remaining formulae to appendix C. The various fits can be acquired in a 
form suitable for inclusion in a Fortran program by contacting one of the authors. 

In Fig. panels A-D, we present the functions Sij(y) for the nuclei 127 I, 129 Xe, 131 Xe, and 
125 Te. The solid lines are for the calculations using the Bonn A Hamiltonian and the dashed 
lines are for the Nijmegen II based Hamiltonian. In order to make comparisons with other 
work easier, we restrict the results to values of y < 2 (q 2 <~ 60000 MeV 2 ). For illustration, 
in Fig. [5] we show the full structure function of 127 I out to y — 10. In Fig. |] panels A-D, 
we show the full structure functions for a pure B [Z\ = 1, Z 2 = Z 3 = Z 4 = 0) for each of 
the nuclei out to a value of y = 2. In these figures each function has been normalized to 
the value S(y = 0) = 1 in order to highlight the similarities and differences in the shapes 
of the structure functions. In our definition of the B we use the older, EMC, values of the 
spin content of the proton. This convention makes it easier to compare our work to previous 
work on 131 Xe ^TJ,Q. The precise values of a p and a n (or a^ and a\) can be found in |22 



the ratio is: a^/ai = 0.297. 

In Figs. H and [5] all of the structure functions have been normalized to S(y = 0) = 1 to 
highlight their similarities and differences. In order to correctly gauge the true differences 
between the various S(q), the different normalizations must be taken into account. This 
is easily done by using eq. (9) and Table |T|. As an example, consider 127 I in panel A of 
Fig. [| and in Fig. |5|. To truly compare the structure functions each of the lines needs to 
be multiplied by a factor such that the ratio at y = is given by Bonn A (Wood-Saxon) : 
Nijmegen II (Wood-Saxon) : Bonn A (Harmonic Oscillator) : Phenomenological (w/ OGM) 
: Single Particle (Harmonic Oscillator) = 1 : 1.42 : 1 : 1/14.7 : 3.47. Similar results can be 
recovered for the other nuclei considered. 

The line labeled Phenomenological above and in the figures requires some explanation. 
This is a shape for a general structure function postulated and used in Ref. []I3"| . It is 
apparent from the figures that this approximation does a reasonable job in reproducing S(q) 
for y < 2. It is clearly inadequate for larger values of the momentum transfer. Below y = 2, 



its shortcomings are also clear but any result derived using this parameterization of S(q) 



for 127 I should not be far off. The overall 127 I axial result of Ref. [13] is another matter, 
as that paper normalizes to the OGM at y = 0. As we have shown, that model severely 
underestimates S(0). 

While the parameterization of S(q) in [[HJ is adequate for I (although we now advocate 
the use of the 127 I structure functions presented in the appendices), it is not applicable for 
all nuclei. The flattening observed in S(q) near y — 1 is the result of higher order multipoles 
becoming important in eq. (8). For 127 I the L = 1,3,5 multipoles all contribute to S(q). 
For small y (<^C 1), the structure function is dominated by the L = 1 multipole. For y > 1, 
all three multipoles contribute and the higher order multipoles dominate. For J = \ nuclei, 
such as 129 Xe and 125 Te, only the L = 1 multipole can contribute. Figs. ||] and § panels B 
and D clearly show that there is no flattening of S(q). Hence, an approximate form like 
that in flTJ] is clearly inappropriate in these cases. In Fig. |6] we show the Bonn A derived 
structure functions for a pure B for all four nuclei. It is obvious that they can not all be fit 
by a single, simple, parameterization. Figs. |] and [5] do show that the pure single particle 
form factor also does an acceptable, but not compelling, job of representing the structure 
functions at all momentum transfers if correctly normalized at y — 0. The correct single 
particle form factor can be easily found by using the tables in the paper by Donnelly and 



Haxton of Ref. ]4C 



Examining the structure functions for 125 Te and 129 Xe in Fig. || illustrates an interesting 
feature. Both of these nuclei are J = \ nuclei with an unpaired neutron. In the ISPSM both 
of these nuclei would be represented by a neutron in the 3s]/2 orbital and have virtually 
identical properties. Table shows that the magnetic moments are quite similar but that 
the distribution of the angular momentum in each nucleus is quite different. This is most 
obvious in the orbital angular momentum L, where the two distributions are quite different. 
Fig. |] reveals that while the structure functions have definite similarities, there are significant 
differences as well. We point all of this out to highlight the fact that seemingly very similar 
nuclei can have very different properties when examined in detail. If precise information on 
the spin distribution of a nucleus is required, detailed calculations must be performed. 

It is also useful to consider differences in S(q) that are the result of different nuclear 
models. S(q) has been calculated for 131 Xe in the context of two other nuclear models, the 
QTDA [B(J and the TFFS [f0J, as well as here. In Fig. [^ we show S(q) for a pure B as a 



function of q 2 for 131 Xe. This figure is meant to be a direct analog of Fig. 2 in Ref. |H| and 



Fig. 3 of Ref. EM. Examining the 3 figures yields some interesting conclusions. All three 



calculations show significant quenching compared to the single particle estimate. The spin 
distribution between the QTDA and TFFS is somewhat different while the full structure 
functions are quite similar. While the values for S ra differ very little between our work 
and the QTDA, the difference in the values of S(0) is almost a factor of 2 between the two 
calculations. Finally, it should be noted that both the QTDA and TFFS calculations of S(q) 
asymptote to the single particle structure function. This is not the case in our calculations, 
which are well below the single particle estimate for all values of q 2 . This can also be seen 
in Ref. |34| where our values of Sij(q) for the Bonn A calculation are compared to those 



of the QTDA calculation. In that comparison, it is apparent that the shell model derived 
structure functions have a much steeper fall off as a function of q 2 . 

Finally, we mention the difference between the structure functions derived using Wood- 
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Saxon wave functions vs. those derived using a harmonic oscillator basis. In Fig. [| panel A 
and Fig. |5| we show the structure functions for 127 I using both sets of basis states. Significant 
differences between the two sets are apparent for extremely high momentum transfers but 
in the range that is most relevant for dark matter detection there is little difference. 

In this section we have discussed the formalism of, and presented our results for, the 
\N axial structure function for several nuclei involved in dark matter detectors. Accurate 
fits which are suitable for use in calculating event rates in detectors are presented in the 
appendices. Several interesting features of the functions have been noted and it is apparent 
that no single simple parameterization of S(q) is suitable for all nuclei. Finally we have 
compared our results to other calculations of 131 Xe structure functions and noted several 
similarities and differences that arise from different nuclear models. 

IV. DISCUSSION 

In this paper we have calculated the full axial response for several heavy nuclei used in 
a number of direct dark matter detection experiments. With this set of structure functions, 
there now exists accurate calculations of the axial xiV response to most, if not all, nuclei 
used as targets in dark matter detectors. We have used the largest model spaces practical 
in conjunction with realistic nuclear Hamiltonians to construct our wave functions. Two 
different nuclear Hamiltonians have been used in order to investigate the sensitivity of our 
results to this particular input. 

The differences in the response due to the two forces is clearly visible in Table [II] and 
Figs. Q-^|. In all cases, reasonable agreement between calculation and experiment for the 
magnetic moment (using free particle ^-factors) is achieved. It is obvious from the table 
that the differences between the two calculations are non-trivial but that they are quite 
a bit smaller than the differences coming from the use of alternate nuclear models. This 
shows that the interaction is not the primary uncertainty in calculations of the \N nuclear 
response. 

We have also attempted to examine the uncertainty due to the nuclear model chosen. 
A number of calculations of 131 Xe's response have been performed. We find that our cal- 
culations are in reasonable agreement with other studies of the spin distribution and finite 
momentum response but distinct differences do exist. In the case of 131 Xe it is not im- 
mediately obvious which calculation is to be preferred. The calculations presented here 
contain more excitations within the model space and use more modern and realistic nuclear 
interactions than the others in the literature. By restricting excitations within this model 
space, the calculations presented in |2{J included excitations out of the space that we worked 



within. Both calculations reproduce the magnetic moment well, with the QTDA calcula- 
tion doing slightly better. (We note that the QTDA model, and a refined version of it, 
have been applied to 127 I and was unable to reproduce the magnetic moment with sufficient 
accuracy JR|.) 

Another improvement incorporated into these calculations of S(q) is the use of Wood- 
Saxon wave functions to evaluate the multipole operators in eq. (8). The Wood-Saxon 
wave functions made a significant difference at extremely high momentum transfers when 
compared to the usual harmonic oscillator wave functions. At the more modest momentum 
transfers typical of "average" neutralinos, the difference is found to be small. 
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Now that these structure functions are available, we hope that they will be used by all 
experiments based upon the materials that we have studied. This will facilitate comparisons 
between different groups. To date, each experiment has used different structure functions 
in their analyses. A first step in this direction has already been taken by one group. The 
most recent analysis of the experiments based in the Gran Sasso laboratory uses the 127 I 
and 129 Xe structure functions presented here [|ll|]. We hope that other groups will follow 
suit so that all future results can be compared on equal footing. 
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TABLES 

TABLE I. The calculated magnetic and quadrupole moments of I compared to experiment 
for calculations using both effective interactions. For the quadrupole moment, effective charges of 
e p = 1.5e and e n = 0.5e have been used. The magnetic moment calculations use the free particle 
^-factors. We also include the ISPSM estimates of the quantities in order to illustrate the quenching 
obtained. 



Observable 


ISPSM 


Bonn A 


Nijmegen II 


Experiment 


Q20 


4.79 
-0.654 


2.775 
-0.577 


3.150 

-0.577 


2.813 
-0.789 
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TABLE II. A comparison of various calculations of the spin distribution of 127 I, 129 Xe, 131 Xe, 
and 125 Te. Bonn A and Nijmegen II are the calculations presented here. OGM is the Odd Group 
Model of @. IBFM is the Interacting Boson Fermion Model of @. TFFS is the Theory of Finite 
Fermi Systems calculation of |g3[- QTDA is the Quasi Tamm-Dancoff Approximation of [£(|. Here, 
as elsewhere in the paper, we take the matrix elements of the operators between the maximally 
stretched states, e.g. (S p ) = (J,Mj = J\S P \J, Mj = J). A blank entry means that the value of 
that particular angular momentum component was not presented in the reference. An entry of 
N/A in the magnetic moment column implies that the experimental magnetic moment was used 
to find the values of spin (S p ) or (S n ) shown. Calculations of the magnetic moment using effective 
g-factors as described in the text are given in parenthesis. 





(S P ) 


(Sn) 


(L P ) 


(Ln) 


1* 


127j 


Experiment 










2.813 


Bonn A 


0.309 


0.075 


1.338 


0.779 


2.775 (2.470) 


Nijmegen II 


0.354 


0.064 


1.418 


0.664 


3.150 (2.790) 


OGM |3C| 


0.07 


0.0 


2.43 


0.0 


N/A 


IBFM |1] 


0.166 


0.01 






N/A 


TFFS H| 


0.15 










ISPSM [32] 


0.5 


0.0 


2.0 


0.0 


4.793 


129 Xe 


Experiment 










-0.778 


Bonn A 


0.028 


0.359 


0.227 


-0.114 


-0.983 (-0.634) 


Nijmegen II 


0.0128 


0.300 


0.372 


-0.185 


-0.701 (-0.379) 


OGM |3C| 


0.0 


0.2 


0.0 


0.3 


N/A 


IBFM |31] 


0.0 


0.2 






N/A 


TFFS H| 




0.25 








ISPSM [3| 


0.0 


0.5 


0.0 


0.0 


-1.913 


131 Xe 


Experiment 










0.69 


Bonn A 


-0.009 


-0.227 


0.165 


1.572 


0.980 (0.637) 


Nijmegen II 


-0.012 


-0.217 


0.215 


1.514 


0.979 (0.347) 


QTDA [20| 


-0.041 


-0.236 


0.026 


1.751 


0.70 


OGM |3C| 


0.0 


-0.18 


0.0 


1.68 


N/A 


IBFM |31] 


0.0 


-0.17 






N/A 


TFFS |3| 




-0.186 








ISPSM H 


0.0 


-0.3 


0.0 


1.8 


1.71 


125 Te 


Experiment 










-0.889 


Bonn A 


0.001 


0.287 


0.077 


0.135 


-1.015 (-0.749) 


Nijmegen II 


-0.0003 


0.323 


0.102 


0.075 


-1.134 (-0.824) 


OGM |3C| 


0.0 


0.23 


0.0 


0.27 


N/A 


IBFM 1 31] 


-0.0004 


0.23 






N/A 


TFFS |3| 




0.22 








ISPSM |32] 


0.0 


0.5 


0.0 


0.0 


-1.913 
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TABLE III. The proton (p) and neutron (n) occupation numbers obtained for each orbital in 
the I calculations. The points to notice are the similarities between the two different interactions 
and the fact that the \hn/ 2 n number is significantly less than 8, the maximum number allowed. 



Orbital 



Bonn A 



Nijmegen II 



I57/2P 

157/2" 



n 



2^5/2 

3s 1/2 p 
3s 1 / 2 n 
2d 3/2 p 



2d 



3/2 



n 



lh 11/2 p 

lh 11/2 n 



1.97979 
7.87440 
0.89648 
5.92205 
0.02859 
1.57985 
0.09511 
2.62799 
0.00004 
5.99571 



1.83023 
7.94902 
0.95545 
5.93768 
0.09023 
1.71861 
0.12374 
2.71644 
0.00034 
5.67825 
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APPENDIX A: 23 Na 

All of the current dark matter detectors which use iodine as a target also use sodium. 
The detectors are large sodium iodide (NaT) crystals ||14| . Since a detailed calculation of the 



axial response of 23 Na has not appeared in the literature, we present one here. The nucleus 
23 Na lies in the middle of the sd shell and therefore is amenable to the same methods applied 
to other srf-shell nuclei. For our calculation we perform the exactly analogous calculation 
to those done for 29 Si in Ref. |22] and 27 A1 in [p5|] ; including the use of harmonic oscillator 



wave functions. The details of the calculation can be found in the above references. 

For 23 Na we use an oscillator parameter of b = 1.6864 fm = (l/117.01)MeV _1 . For our 
adopted maximum halo velocity of v = 700 km/sec we have y ma x = 0.1875. A breakdown of 
the angular momentum along with a comparison of the measured and calculated magnetic 
moments is presented in Table |V|; agreement is excellent. Table |IV| also shows a significant 



difference in S p from that predicted in the OGM. Finally, in the following equations we 
present fits to the structure functions Sij(q) as 3rd order polynomials in y which are highly 
accurate to values well past y ma x- 

TABLE IV. The decomposition of the angular momentum for 23 Na along with the calculated 
and experimental magnetic moments. 





(S P > 


(Sn) 


<L P > 


(Ln> 


A* 


23 Na 


Experiment 

Calculation 

OGM 


0.2477 
0.1566 


0.0198 
0.0 


0.9117 
1.3434 


0.3206 
0.0 


2.218 
2.2196 

N/A 



S 00 (y) = 0.0379935 - 0.174341?/ + 0.378299?/ 2 - 0.342962?/ 3 (Al) 

S 01 (y) = 0.0646525 - 0.350289?/ + 0.910031?/ 2 - 0.985833?/ 3 (A2) 

S n (y) = 0.0275013 - 0.169641y + 0.507868?/ 2 - 0.617985y 3 (A3) 
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APPENDIX B: THE ABBREVIATED STRUCTURE FUNCTIONS 

The fits to S(q) in this appendix are only valid for y < 1. The fits are presented as tables 
of the coefficients of 6th order polynomials in y: SV,(g) = J2t=o CkU k - The first column gives 
the order of y k then next 3 columns give the corresponding values of the C& for Soo, S i, an d 
Sn for the Bonn A calculation. The last 3 columns present the results for the Nijmegen II 
calculation in the same manner. 









TABLE V. 


127j 










Bonn A 






Nijmegen II 






5*00 


5*01 


Sn 


Soo 


5*01 


Su 


1 


0.0982724 


0.119851 


0.0365375 


0.116548 


0.161931 


0.0562404 


y 


-0.675013 


-0.843567 


-0.262676 


-0.792274 


-1.14026 


-0.408512 


y 2 


2.13531 


2.73535 


0.875115 


2.49846 


3.71441 


1.37775 


y 3 


-3.7595 


-4.93029 


-1.61455 


-4.38312 


-6.71583 


-2.57019 


y 1 


3.77735 


5.05806 


1.69076 


4.38495 


6.89384 


2.70866 


y 5 


-2.0091 


-2.73609 


-0.930164 


-2.32223 


-3.72586 


-1.4945 


y 6 


0.435566 


0.60084 


0.206944 


0.501504 


0.817068 


0.332885 








TABLE VI. 


129 Xe 










Bonn A 






Nijmegen II 






•Soo 


Soi 


Su 


Soo 


Soi 


Su 


l 


0.0712796 


-0.121583 


0.0518388 


0.0464592 


-0.0853234 


0.0391694 


y 


-0.480418 


0.874546 


-0.394855 


-0.313776 


0.614961 


-0.299123 


y 2 


1.47263 


-2.83165 


1.34334 


0.965631 


-1.98471 


1.00873 


y 3 


-2.53226 


5.09221 


-2.51522 


-1.6666 


3.54959 


-1.86483 


v 1 


2.49681 


-5.19757 


2.64796 


1.64774 


-3.60225 


1.93996 


yb 


-1.30712 


2.79235 


-1.4557 


-0.864196 


1.92573 


-1.05644 


y e 


0.279589 


-0.60881 


0.322793 


0.185069 


-0.418234 


0.232626 








TABLE VII. 


131 Xe 










Bonn A 






Nijmegen II 






5*00 


Soi 


Su 


Soo 


5*01 


Sn 


l 


0.0295866 


-0.0544505 


0.0250499 


0.0277038 


-0.0497326 


0.0223178 


y 


-0.185155 


0.36762 


-0.181162 


-0.175382 


0.338942 


-0.162659 


y 2 


0.593387 


-1.18133 


0.593168 


0.560377 


-1.10015 


0.542687 


y 3 


-1.03518 


2.05291 


-1.03886 


-0.996936 


1.97087 


-0.98921 


2/ 4 


1.00492 


-1.98269 


1.00706 


1.01 


-1.99963 


1.01495 


y'' 


-0.507773 


0.996715 


-0.50709 


-0.540224 


1.06809 


-0.54588 


y 6 


0.103658 


-0.202596 


0.103134 


0.11739 


-0.231591 


0.118858 
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TABLE VIII. 


125 Te 










Bonn A 






Nijmegen II 






5*00 


Sbi 


Sn 


Soo 


Soi 


Sn 


1 


0.0396831 


-0.0788638 


0.0391772 


0.049567 


-0.0993273 


0.0497519 


y 


-0.271174 


0.572717 


-0.30043 


-0.342464 


0.731525 


-0.387508 


y 2 


0.869383 


-1.90069 


1.03775 


1.06657 


-2.39301 


1.32901 


y 3 


-1.56951 


3.46977 


-1.94604 


-1.85469 


4.32285 


-2.49519 


y 4 


1.61835 


-3.5546 


2.02635 


1.84644 


-4.42823 


2.6334 


y 5 


-0.879731 


1.90199 


-1.09438 


-0.976196 


2.39046 


-1.45353 


y 6 


0.194048 


-0.411614 


0.238043 


0.210989 


-0.524466 


0.32414 



APPENDIX C: THE FULL STRUCTURE FUNCTIONS 



The fits to S(q) in this appendix are good for all values of y < 10. The fits are presented 
as tables of the coefficients of 8th order polynomials in y plus a term included to mimic 
the Goldberger-Trieman term present in the longitudinal multipole p0 |, all multiplied by a 
factor of exp(— 2y) : Sij(q) = (J2k=o^ky k + Cgj^ r )e~ 2y . The first column gives the order of 



y then next 3 columns give the corresponding values of the Ck for Soo, Soi, and Sn for the 
Bonn A calculation. The last 3 columns present the results for the Nijmegen II calculation 
in the same manner. An example of the table's use can be found by comparing Eq. (11) to 
the entries for Sqq in the Bonn A portion of the 127 I table. 









TABLE IX. 


127j 










Bonn A 






Nijmegen II 




x{e~ 2 y) 


•Soo 


Soi 


Sn 


Soo 


Soi 


S\i 


1 


0.0983393 


0.11994 


0.0365709 


0.11663 


0.162054 


0.056287 


y 


-0.489096 


-0.618424 


-0.194994 


-0.572149 


-0.836288 


-0.303825 


y 2 


1.1402 


1.50893 


0.504876 


1.33797 


2.05944 


0.794783 


y 3 


-1.47168 


-2.07367 


-0.747451 


-1.72517 


-2.83193 


-1.17027 


</ 


1.1717 


1.77307 


0.704334 


1.37742 


2.39726 


1.06373 


if 


-0.564574 


-0.903597 


-0.393018 


-0.669986 


-1.21214 


-0.571342 


if 


0.158287 


0.26002 


0.121881 


0.190522 


0.348612 


0.172197 


y 7 


-0.0238874 


-0.0387025 


-0.0191881 


-0.0291803 


-0.0521813 


-0.0266165 


if 


0.00154252 


0.00235675 


0.00121021 


0.0019081 


0.00320731 


0.00166238 


i 

l+j/ 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



TABLE X. 129 Xe 







Bonn A 






Nijmegen II 




x{e~ 2 y) 


Soo 


Soi 


Su 


Soo 


Soi 


Su 


1 


0.0713238 


-0.12166 


-2.05825 


0.046489 


-0.0853786 


-1.28214 


y 


-0.344779 


0.644351 


1.80756 


-0.225507 


0.453434 


1.09276 


y 2 


0.755895 


-1.52732 


-1.27746 


0.499045 


-1.06546 


-0.712949 


if 


-0.933448 


2.02061 


0.654589 


-0.622439 


1.3867 


0.314894 


y 4 


0.690061 


-1.57689 


-0.221971 


0.46361 


-1.0594 


-0.0835104 


y 5 


-0.302476 


0.723976 


0.0454635 


-0.20375 


0.47576 


0.0105933 


if 


0.0765282 


-0.190399 


-0.00425694 


0.0510851 


-0.122077 


0.000233709 


if 


-0.0103169 


0.0263823 


-0.000136779 


-0.00670516 


0.0164292 


-0.000243292 


if 


0.000573919 


-0.00148593 


0.00004396 


0.00035659 


-0.000894498 


0.0000221666 


i 

l+y 


0.0 


0.0 


2.11016 


0.0 


0.0 


1.32136 



TABLE XL 131 Xe 







Bonn A 






Nijmegen II 




x{e~ 2 y) 


Soo 


Sbi 


Su 


Soo 


Soi 


Su 


1 


0.0296421 


-0.0545474 


0.0250994 


0.0277344 


-0.0497844 


0.0223447 


y 


-0.133427 


0.271757 


-0.137716 


-0.124487 


0.247247 


-0.122063 


y 2 


0.377987 


-0.723023 


0.366609 


0.328287 


-0.632306 


0.319493 


if 


-0.579614 


1.0545 


-0.53851 


-0.481399 


0.896416 


-0.466949 


y 4 


0.578896 


-0.971333 


0.492545 


0.475646 


-0.816445 


0.428767 


y 5 


-0.345562 


0.538422 


-0.269903 


-0.285177 


0.452352 


-0.236789 


if 


0.115952 


-0.168988 


0.0836943 


0.0968193 


-0.142686 


0.0740837 


if 


-0.0201178 


0.027416 


-0.0133959 


-0.0170957 


0.0233463 


-0.0119668 


f 


0.00141793 


-0.00180527 


0.000868668 


0.00123738 


-0.00156293 


0.000787042 


i 

1+2/ 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 
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TABLE XII. 125 Te 







Bonn A 






Nijmegen II 




x{e~ 2 y) 


Soo 


Sbi 


Sn 


5*00 


Soi 


Si i 


1 


0.0397091 


-0.0789431 


0.0392236 


0.0495946 


-0.0993873 


-1.92941 


y 


-0.196101 


0.42738 


-0.229376 


-0.247766 


0.54303 


1.68075 


y 2 


0.472653 


-1.09331 


0.622146 


0.547656 


-1.28816 


-1.16336 


y ?> 


-0.650229 


1.55324 


-0.922531 


-0.665532 


1.67206 


0.586501 


y 4 


0.541926 


-1.28933 


0.784648 


0.474621 


-1.26883 


-0.207302 


y 5 


-0.264563 


0.618441 


-0.382445 


-0.199442 


0.56728 


0.0514094 


if 


0.074891 


-0.16964 


0.105709 


0.0481866 


-0.145438 


-0.00869728 


y 7 


-0.0114632 


0.0248165 


-0.0154157 


-0.00616326 


0.0195887 


0.000870366 


j/ 


0.000749022 


-0.00152108 


0.000928651 


0.000322728 


-0.00106519 


0.0000354095 


i 

l+y 


0.0 


0.0 


0.0 


0.0 


0.0 


1.97923 
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FIGURES 

FIG. 1. A visual description of the 127 I model space. See the text for specific details of the 
construction of the space. The other nuclei studied use the same jump assignments and SPEs. 
The left 1*77/2 SPE is that used for the Bonn A interaction and the one on the right is used by the 
Nijmegen II interaction. 

FIG. 2. The three spin structure functions Sij(q) where i,j = 0, 1 for the 4 nuclei considered. 
The results using both effective interactions are plotted. Accurate fits to these structure functions 
can be found in Appendices B and C. A) Sij(q) for 127 I; the ordering is: 5rji > 5oo > 5n f° r each 
force. B) Sij(q) for 129 Xe; the ordering is: 5oo > 5n > 5oi f° r each interaction. C) Sij(q) for 
131 Xe; the ordering is: Soo > S11 > 5rji f° r each interaction. D) Sij(q) for 125 Te; the ordering is: 
Soo > "Sii > S01 ^ or each interaction. 

FIG. 3. Another view of panel A in Fig. |2| Here we have extended Sij(q) out to y = 10 
and chopped off much of the initial fall off from 5^(0) in order to highlight the similarities and 
differences between the two sets of structure functions. 

FIG. 4. The spin structure function S(q) for a pure B (ao/ a i = 0.297) for the 4 nuclei con- 
sidered. Wood-Saxon (W.S.) wave functions have been used. The results using both effective 
interactions are plotted. Additionally, the pure single particle estimate of S{q) with harmonic 
oscillator (H.O.) wave functions is included for comparison. All structure functions have been 
normalized to 5(0) = 1 in order to better compare their intrinsic shapes. To truly compare the 
differences the functions need to be normalized using Eq. (9) and the values in Table O. A) S(q) 
for T. Also included for comparison is the S(q) used in |l!| and the results for the Bonn A 
interaction using H.O. wave functions. B) S{q) for 129 Xe. C) 5(g) for 131 Xe. C) 5(g) for 125 Te. 

FIG. 5. Another view of panel A in Fig. ||. Here we have extended 5(g) out to y = 10 and 
chopped off much of the initial fall off from 5(0) in order to highlight the similarities and differences 
between the various structure functions. 

FIG. 6. The Bonn A calculations of S(q) for a pure B for all 4 nuclei compared. Note the very 
large differences near y = 1 between the nuclei with J = ^ ( 129 Xe and 125 Te) and those with larger 
J. 

FIG. 7. The 131 Xe structure function for a pure B. The single particle structure function has 
been normalized to 5(0) = 1. The Bonn A and Nijmegen II calculations have been correctly 
normalized relative to the single particle model. This figure is a direct analog of, and should 
be compared to, Fig. 2 of Ref. [EQ] and Fig. 3 of Ref. |^4j. The major differences between 
these calculations and those of p3| and [44] are: We find 5(0) ~ 0.4-0.5 vs. their values of 



5(0) ~ 0.25 and both of the other model's structure functions asymptote to the single particle 
model for q 2 > 0.02 GeV 2 while these calculations stay well below the single particle model. 
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